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Abstract The failure behavior of wood with a short crack 
was examined by conducting the single-edge-notched 
bending tests of a radial-longitudinal system on Agathis 
specimens. In the test, the mode I critical stress intensity 
factor was measured, and its validity was checked by the 
result from double cantilever beam testing method. The 
mode I critical stress intensity factor decreased when 
the crack length approached zero. With crack length 
correction, a constant critical stress intensity factor was 
obtained over a wide range of crack length including crack- 
free specimen. 
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Introduction 

In actual cases as timber engineering and wood drying, the 
crack lengths in wood are often short. In a previous study 
[1], the single-edge-notched bending (SENB) tests of tan¬ 
gential-radial (TR) system were conducted using speci¬ 
mens of Agathis (Agathis sp.) with varying the crack 
lengths that were smaller than those adopted in several 
existing studies involving the single-edge-notched tension 
or bending tests [2-8]. When considering an additional 
crack length ahead of the crack tip, the mode I critical 
stress intensity factor was effectively predicted indepen¬ 
dently from the initial crack length. In addition, the 
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nominal bending strength could also be predicted in both 
cases of with and without a crack based on the concept for 
the additional crack length. 

In practical operations such as timber engineering and 
drying, a crack induced in solid wood often propagates 
along the fiber, so it is important to investigate whether the 
concept described above is valid for the system including 
the fiber direction. In this study, the radial-longitudinal 
(RL) system SENB tests were conducted, and the validity 
of the results was examined by comparing test data with 
those obtained from double cantilever beam (DCB) tests. 


Theory 

Three-point SENB test analyses 

Figure 1 shows a schematic diagram of the three-point 
SENB test. The specimen which had a crack of length a at 
its center was supported with a span of S , and the load P 
was applied at the mid-span. The initial crack length is 
defined as a. In this loading condition, the nominal bending 
stress <r n is derived from the elementary beam theory as 
follows: 

_ 3 SP 

a " 2BW 2 ’ ' 

where B and W are the beam width and depth, respectively. 
This notation is applicable to crack-free specimens. When 
a crack-free specimen is bent, the failure-by-bending 
moment is induced when o n reaches its critical value cr nc , 
which is usually defined as the bending strength of the 
material. 

The fracture behavior of a material with a long crack 
should be analyzed using the mode I stress intensity factor 
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Fig. 2 Schematic diagram of double cantilever bending (DCB) test 


Fig. 1 Schematic diagram of the three-point single-edge-notched 
bending (SENB) test 
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Ki or the energy release rate G h each of which is derived 
from fracture mechanics theory. In SENB testing of spec¬ 
imens with a short crack, however, it is difficult to measure 
Gi appropriately, because the loading-line compliance does 
not vary with the crack length [1]. In contrast, Kj is more 
easily measured using an approximating equation. 

The mode I stress intensity factor Kj is derived using the 
following equation: 

K\ = <x nX /ia/Q^, (2) 


where f(a/W) is the crack geometry factor. The fracture is 
initiated when Ki reaches the critical stress intensity factor, 
defined as AT Ic , which is obtained by substituting a nc for a n 
in Eq. (2). Gross and Srawley derived a f(aJW)-aJW 
relationship for isotropic material as follows [9]: 



= 1.09- 1.735 f—3 

\WJ 

+ 8 - 20 (w) 2 - 1418 (w ) i+1457 




There was a concern that the orthotropy of the wood 
might influence the crack geometry factor. However, 
Eq. (3) was appropriate as the crack geometry factor for 
the TR system SENB specimens with the configuration 
similar to this study according to the virtual crack closure 
technique (VCCT) applied to the results obtained from the 
finite element analyses [1]. Moreover, the validity of 
Eq. (3) was verified for wood with a strong orthotropy, 
such as TL system fracture mechanics testing analyses 
[10]. Hereafter, the mode I stress intensity factor is 
calculated using the crack geometry factor represented by 
Eq. (3). 


where C L is the compliance obtained from the load- 
deflection relationship, and C s is the compliance obtained 
from the load-longitudinal strain relationship, which was 
measured at / = a/2. 

Here, the energy release rate was converted to the stress 
intensity factors as follows: 

K, = ^ (5) 

where 


c i 



Ex 1 (Ex_ 
Ey 2 \G X} ; 




in which E x and E y are Young’s moduli in the x and 
y directions, respectively, and G xy and v*y are the shear 
modulus and Poisson’s ratio in the xy plane, respectively 
[12]. In this study, xy plane denotes RL plane. The elastic 
modulus values were obtained by compression tests and 
vibration tests [1], and are listed in Table 1. 


Experimental 

Materials 

Agathis (Agathis sp.) lumber, with a density of 
480 ±10 kg/m 3 at 12 % moisture content (MC), was used 
for the tests. When examining the fracture mechanics 
properties of the radial-longitudinal system, the influence 
of the annual rings is often significant [8]. Agathis is a 
tropical wood species with no significant growth rings; 


Double cantilever beam test analyses 

Double cantilever beam test in the RL system (Fig. 2) was 
conducted to examine the validity of SENB testing results. 
According to beam theory and fracture mechanics, the 
mode I energy release rate can be obtained as follows [11]: 


Table 1 Elastic constants used for the conversion to stress intensity 
factor 


E x (GPa) 

E y (GPa) 

G xy (GPa) 


14.6 

0.87 

1.0 

0.44 


x and y directions correspond to the longitudinal and radial directions 
of wood, respectively 


<£) Springer 















































J Wood Sci (2012) 58:453-458 


455 


hence, the influence of annual rings could be ignored in the 
analysis. The lumber had no defects such as knots or grain 
distortions, so that the specimens cut from it could be 
regarded as small and clear. The lumber was stored for 
several months in a room at a constant temperature of 
20 °C and a relative humidity of 65 % before the test and 
was confirmed to be in an air-dried condition. These con¬ 
ditions were maintained throughout the tests. The equilib¬ 
rium MC was approximately 12 %. 

Single-edge-notched bending tests 

All of the specimens were cut from the lumber described 
above so that they were long-matched to the dimensions 
shown in Table 2. The numbers of specimens were 50,25, and 
15 for specimen A, B, and C, respectively. For the cracked 
specimens, a crack was produced along the radial direction in 
the RL system. The crack was first cut with a band saw 
(thickness = 1 mm), it was then extended ahead of the crack 
tip using a razor blade. The specimen was supported by the 
span shown in Table 2, and a load was applied to the specimen 
at a crosshead speed of 1 mm/min for the test with span lengths 
of 60 and 120 mm, a crosshead speed of 2 mm/min was used 
for the test with a span length of 240 mm. The test was con¬ 
ducted until the load markedly decreased. 

Figure 3 shows a typical load-loading-line deflection 
relation obtained by the SENB test. The load increases lin¬ 
early until it is maximum, and then drops immediately 
without displaying a nonlinear phase. This trend is similar to 
that obtained from the SENB tests of the specimen with the 
TR system [1]. In this research, the critical load for crack 
propagation P c was determined as the maximum load, as the 
same manner with that in the previous study [1]. By substi¬ 
tuting P c into Eq. (1), the nominal bending strength a nc 
corresponding to initial crack length a was obtained [1]. 
Then, the mode I critical stress intensity factor Ki c corre¬ 
sponding to the initial crack length was obtained by substi¬ 
tuting er nc and f(a/W) represented as Eq. (3) into Eq. (2). 

Double cantilever beam test 

All of the specimens were cut from the lumber mentioned above 
so that they were side matched into the initial dimensions of 
30 mm (radial direction) x 15 mm (tangential direction) x 
315 mm (longitudinal direction). The number of specimens was 


Table 2 Specimen configurations used for the SENB tests 


Specimen 

type 

Span 

S (mm) 

Depth 

W (mm) 

Width 

B (mm) 

Crack length 
a (mm) 

A 

60 

15 

7.5 

0, 0.5, 1, 2, 4 

B 

120 

30 

15 

0, 0.5, 1, 2, 4 

C 

240 

60 

30 

0, 0.5, 1, 2, 4 



Fig. 3 Typical example of the load-loading-line deflection relation 
by the SENB test 

six. The crack was first cut along the longitudinal direction in the 
longitudinal-radial plane with a band saw (thickness = 1 mm), 
and then it was extended 1 mm ahead of the crack tip using a 
razor blade so that the crack length would be as that mentioned 
below. After cutting the crack, the loading blocks of MDF with 
the dimensions of 30 mm in length, 30 mm in height, and 
15 mm in thickness were bonded by epoxy resin on the upper 
and lower cantilever portions symmetrically as shown in Fig. 2. 
Crack length, which was defined as the distance from the line of 
load application to the crack tip, was 150 mm. Load P was 
applied to the specimen by pins through universal joints at a 
crosshead speed of 5 mm/min until the load markedly 
decreased. The total testing time was about 10 min. 

The loading-line displacement 3 was measured by the 
crosshead travel since it was confirmed that the machine 
compliance was small enough to be ignored [7], whereas the 
longitudinal strain s was measured by a strain gauge (gauge 
length = 2 mm; FLA-2-11, Tokyo Sokki Co., Tokyo) bon¬ 
ded at the midpoint between the loading line and crack tip on 
the top cantilever portion (/ = all). The C L and C s values 
were obtained from the linear portions of P-3 and P-s 
relationships, respectively [11]. The critical load for crack 
propagation P c was defined as that at the onset of nonlinearity 
in the load-loading-line displacement curve [7, 13, 14]. 

Results and discussion 

Three-point SENB test 

Figure 4 shows the relationship between the nominal bending 
strength er nc and initial crack length a. According to the results 
for metals [15,16], the nominal bending strength of a cracked 
specimen will approach the bending strength of a crack-free 
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Fig. 4 Relationship between the nominal bending strength a nc and 
crack length a by the SENB test 


specimen when the crack length is decreased. From the results 
obtained here, however, the value of cr nc is markedly smaller 
than that of the crack-free specimen, and it decreases as the 
crack length increases. This tendency is similar to that for the 
results obtained using the specimens with the TR system [1], 
and it suggests that fracture mechanics theory is essential for 
analyzing the failure behavior of solid wood with a crack even 
when the crack length is short. 

Figure 5 shows the relationship between the mode I 
critical stress intensity factor 7f Ic and the initial crack 
length a. The K ic value obtained by the DCB test, the 
details of which are mentioned below, is also shown in 
Fig. 5. The AT Ic value can be regarded to be a material 
parameter independently determined from the initial crack 
length. Nevertheless, the AT Ic value is not constant and it 
significantly increases as the initial crack length increases. 

The dependence of 7f Ic value on the initial crack length a 
can be reduced by introducing an additional crack length A 
into Eq. (2) [1, 17]: 

Ki = <T n a/t i(a + A)/ = OnVna'f (^) , (7) 

where a! is defined as the corrected crack length. Table 3 
shows the value of the additional crack length A 



0 1 2 3 4 5 

Crack length a (mm) 


corresponding to each specimen type. The value of A was 
determined as follows: (1) the values of K ic were calculated 
under various values of A in Eq. (7), (2) the standard 
deviations of K ic corresponding to each crack length were 
obtained, (3) the sums of the standard deviations, defined as 
s D , corresponding to each A were obtained and compared. 
The value of A in Table 3 was determined as that from 
which the smallest value of s D was derived. Figure 6 shows 
the relationship between the mode I critical stress intensity 
factor Ki c and the corrected crack length a !. By conducting 
the crack length correction, the dependence of Ki c on the 
crack length is effectively reduced. The concept for intro¬ 
ducing the additional crack length is also applicable for the 
crack-free specimen. Although the dimensions of the 
specimens vary significantly, the variation of the additional 
crack length values is relatively small. These additional 
crack length values are similar to those obtained from the 
tests of the specimens with the TR system [1]. In a previous 
study, we examined the mode II fracture behavior of Aga- 
this by asymmetric four-point bending tests of a specimen 
with a short crack in the TL system, which were cut from 
the same lumber used in this investigation [17]. The test 
results were analyzed with the same manner conducted 
here, and we found that the additional crack length values 
obtained in this study are shorter than that obtained from the 
mode II tests [17]. Table 4 shows the average values of Ki c 
obtained with and without correcting the crack length. 

The relationships between a nc and a can be conven¬ 
tionally predicted by substituting the values of 7f Ic without 


Table 3 Additional crack length corresponding 
type 

to each specimen 

Specimen 

Additional crack 

type 

length A (mm) 

A 

0.25 

B 

0.33 

C 

0.33 



Fig. 5 Relationship between the mode I critical stress intensity factor 
Ki c and crack length a 


Fig. 6 Relationship between the mode I critical stress intensity factor 
Kic and corrected crack length a' 
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Table 4 Mode I critical stress intensity factor Ki c obtained with and 
without correcting the crack length 


Specimen type 

K lc (MPa m 1/2 ) 


Crack length uncorrected 

Crack length corrected 

A 

0.35 ± 0.052(a) 

0.38 ± 0.037(b) 

B 

0.36 ± 0.069(a) 

0.39 ± 0.051(b) 

C 

0.37 db 0.058(a) 

0.41 ± 0.047(b) 


Results are averages ± SD. The K ic values listed are obtained by 
averaging the values corresponding to the initial crack length a of 0.5, 
1, 2, 3, and 4 mm in (a), and a = 0, 0.5, 1, 2, 3, and 4 mm in (b) 



Fig. 7 Comparison of nominal bending strengths experimentally 
obtained and predicted on the basis of fracture mechanics theory 

the crack length correction in Table 4 and f(a/W) into 
Eq. (2) [15]. In addition, the (J nc -a relationships are also 
predicted by substituting the values of 7f Ic with the crack 
length correction in Table 4, A in Table 3, and f((a + A)/ 
W) into Eq. (7). This method is our proposal. Figure 7 
compares the predicted and experimentally obtained 
cT nc -a relationships. When the crack length is not corrected, 
the predicted strength increases markedly when the crack 
length approaches zero. When the crack length is corrected, 
however, the (J nc -a relationship is effectively predicted 
over the full range of crack lengths, including when a = 0. 



Fig. 8 A relationship between strength and crack length for wood 

This is the same result as those in the previous works [1, 17]. 
Therefore, the failure behavior is probably described by 
fracture mechanics theory even when the specimen has no 
crack as shown in Fig. 8 in comparison with the failure 
behavior of metal [15, 16]. 

For wood as quasi-brittle material, so-called fracture 
process zone influence [18-24] might not be so significant 
because the nonlinearity in the load-loading-line deflection 
relation is not significant as Fig. 3 shows. Natural cracks 
inherently contained in wood might be significant to the 
phenomenon obtained in this study. Further researches 
including the microscopic observation will reveal this 
phenomenon. 

Comparison with DCB test results 

The average Ki c value by DCB test was 0.42 db 0.02 
MPa m 1/2 . The average K ic values by SENB are 0.35 db 
0.058 MPa m " for 0.5-4 mm crack length specimens 
before correction and 0.39 db 0.044 MPa m 1/2 for all 
specimens including crack-free specimens after additional 
crack length correction. The results obtained from the 
SENB tests approached to DCB results after the correction. 
By the t test between the results obtained from the DCB 
and SENB tests, the probability value (P value) was 0.011 
before correction and the difference was significant at 5 % 
level. In contrast, the P value was 0.19 after correction, so 
the difference was not significant at 5 % level. Although 
the DCB specimens used here were not longitudinally 
matched to the SENB specimens but were only in the same 
lumber, the results may show the validity of introducing 
the additional crack length. 

There have been several approaches on the application 
of fracture mechanics to the crack free specimens [25, 26]. 
Further researches are required for revealing the failure 
behavior of solid wood based on our proposal as well as the 
examinations described above. 
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